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ABSTRACT 

In this work high spatial resolution mid-infrared (MIR) 12 fim continuum imaging of low-luminosity active galactic nuclei (LLAGN) 
obtained by VLT/VISIR is presented. The goal of this investigation is to determine if the nuclear MIR emission of LLAGN is consistent 
with the existence of a dusty obscuring torus, the key component of the unification model for AGN. Based on available hard X-ray 
luminosities and the previously-known tight correlation between the hard X-ray and 12yum luminosities, a sample of 17 nearby 
LLAGN without available VISIR N-band photometry was selected. Combined with archival VISIR data of 9 additional LLAGN 
with available X-ray measurements, the dataset represents the bulk of southern LLAGN currently detectable from the ground in the 
MIR. Of the 17 observed LLAGN, 7 are detected, while upper limits are derived for the 10 non-detections. Thus, the total number 
of AGN detected with VLT/VISIR increases to more than 50. All detections except NGC 3125 appear point-like on a spatial scale of 
~ 0.35". The detections do not significantly deviate from the known MIR-X-ray correlation but extend it by a factor of ~ 10 down to 
luminosities < 10 41 erg/s with a narrow scatter (tr = 0.35 dex, Spearman Rank p = 0.92 ). The latter is dominated by the uncertainties 
in the X-ray luminosity. Interestingly, a similar correlation with comparable slope but with a normalization differing by ~ 2.6 orders 
of magnitude has been found for local starburst galaxies. In addition, the VISIR data are compared with lower spatial resolution data 
from Spitzer/IRS and IRAS. By using a scaled starburst template SED and the PAH 1 1.3 /im emission line the maximum nuclear star 
formation contamination to the VISIR photometry is restricted to < 30% for 75% of the LLAGN. Exceptions are NGC 1097 and NGC 
1566, which may possess unresolved strong PAH emission. Furthermore, within the uncertainties the MIR-X-ray luminosity ratio is 
unchanged over more than 4 orders of magnitude in accretion rate. These results are consistent with the existence of the dusty torus 
in all observed LLAGN, although the jet or accretion disk as origin of the MIR emission cannot be excluded. Finally, the fact that the 
MIR-X-ray correlation holds for all LLAGN and Seyferts makes it a very useful empirical tool for converting between the MIR and 
X-ray powers of these nuclei. 

Key words, galaxies: active - galaxies: nuclei - infrared: galaxies - X-rays: galaxies - accretion, accretion disks 



1. Introduction 

Low-luminosity active galactic nuclei (LLAGN; L2-iokeV 55 
10 42 erg/s) have drawn increasing attention in recent years for 
various reasons: it has been realized that this population is not 
only the most common of all active galaxies (AGN) - about two- 
thirds of all AGN are LLAGN (see Ho 2008) - but they may 
also represent AGN in the initial or final stages of their active 
phases. The underlying idea is that every galaxy harbors a black 
hole (BH) in its center that grows and evolves to a supermassive 
BH during many phases of nuclear activity. The key components 
describing AGN in the so called unification scheme (Antonucci 
1993) such as the accretion disk and the dusty obscuring torus 
very likely vanish during non-active phases and are built-up 
again when new material reaches the galactic nucleus due to any 
kind of larger scale disturbance. But it is unknown how the pro- 
cesses of formation and evolution of both are connected. This 
is partially because the physical nature and morphology of the 
disk and torus is still not exactly known. However, high spa- 
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tial resolution mid-infrared (MIR) observations made progress 
in the modeling of the torus possible. They favor a clumpy dis- 
tribution of optically thick clouds with a small filling factor (e.g. 
Nenkova et al. 2002; Dullemond & van Bemmel 2005; Honig 
et al. 2006; Schartmann et al. 2008; Nenkova et al. 2008; Honig 
& Kishimoto 2010). Apart from the fact that the mechanism, 
which maintains the geometrical thickness of the torus, is un- 
known, there are various theoretical works predicting its dis- 
appearance below certain bolometric luminosities or accretion 
rates (Elitzur & Shlosman 2006; Honig & Beckert 2007). But 
observational proof of this disappearance is still lacking. 

The putative dusty torus of AGN is best observed in the MIR, 
where the thermal emission from the warm dust is strongest. 
Ground-based MIR instruments like VLT/VISIR (Lagage et al. 
2004) are best suited to perform these observations due to their 
high spatial resolution. And an increasing number of AGN stud- 
ies demonstrate their power to isolate the AGN from surround- 
ing starbursts (e.g. Gorjian et al. 2004; Galliano et al. 2005; 
Haas et al. 2007; Siebenmorgen et al. 2008; Ramos Almeida 
et al. 2009). Furthermore, mid-infrared interferometric observa- 
tions were able to resolve the dusty torus in nearby AGN and 
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showed that it is only a few parsec in size (e.g. Jaffe et al. 2004; 
Raban et al. 2009; Tristram et al. 2009; Burtscher et al. 2009) 
and scales by the square root of the AGN luminosity (Tristram 
& Schartmann 2011). 

On the other hand, past low-resolution MIR observations led 
to the discovery of a correlation between the MIR and X-ray 
luminosities (Krabbe et al. 2001; Lutz et al. 2004). Here, the 
observed nuclear MIR emission is believed to represent the re- 
processed UV/X-ray radiation from the dusty clouds inside the 
torus. This conclusion is strongly supported by fitting of high- 
spatial resolution spectrophotometry with clumpy torus models 
(Honig et al. 2010; Alonso-Herrero et al. 2011). 

Then, recently gathered VISIR photometry enabled the best 
estimate of AGN emission around 12 micron, free of stellar 
contamination, resulting in a strong correlation with the intrin- 
sic X-ray 2-10keV luminosity (Horst et al. 2006, 2008, 2009; 
Gandhi et al. 2009), which is well-described by a power-law 
(Lmir k L^) with a slope b ~ 1 . Similar results have also been 
found with other ground-based high spatial resolution MIR in- 
struments (e.g. Levenson et al. 2009). The luminosity correlation 
for the different AGN types is indistinguishable within the mea- 
sured precision, and it can be used to e.g. assess the bolometric 
IR luminosity of AGN when no high resolution IR data is avail- 
able (e.g. Mullaney et al. 201 1). 

However, the lowest luminosity objects of the observed sam- 
ples (in e.g. Gandhi et al. 2009, hereafter G+09), in particular 
the 3 LINERs, show a tendency to deviate from this correlation: 
all objects with an X-ray luminosity L2-iokeV < 10 42 erg/s ex- 
hibit a MIR-excess (or an X-ray deficit) of ~ 0.3 dex. The total 
number of included LLAGN was too small, to give statistical 
significant evidence for a real deviation. While such a deviation 
could indicate a change in the physical structure of the AGN at 
low luminosities, one first has to exclude other causes as e.g. sig- 
nificant MIR emission from dust heated by circum-nuclear star 
formation regions. In particular at low luminosities, the latter 
can be of comparable MIR brightness with respect to the AGN 
or even brighter. Furthermore, the proposed collapse of the dusty 
torus (see above) would lead to a deficit of MIR emission, oppo- 
site to the indicated deviation. Instead, a separate AGN compo- 
nent would have to explain the MIR-excess (or X-ray-deficit) at 
low luminosities, if real. One candidate would be the jet, which 
is known to emit copious amounts of synchrotron emission over 
the whole observable wavelength range. In particular in LINERs 
(short for low-ionization nuclear emission-line region, a sub- 
class of LLAGN), the jet is assumed to be much more power- 
ful than in Seyferts, explaining relative radio-brightness of these 
optically faint objects. For example, indications of significant 
jet contribution to the MIR emission of LLAGN was found in 
Hardcastle et al. (2009). In fact for NGC 4486, a well-known 
LINER, the observed MIR emission can be explained only with 
synchrotron emission of different jet components (Perlman et al. 
2007). 

However, the inner accretion disk might also emit MIR 
synchrotron emission, when in the so-called radiatively ineffi- 
cient accretion flow (RIAF) or advection dominated accretion 
flow (ADAF) mode. The latter was theoretically predicted by 
Narayan & Yi (1994) to occur at low accretion rates: the stan- 
dard geometrically thin, optically thick disk (Shakura & Sunyaev 
1973) transforms into a geometrically thick, optically thin accre- 
tion disk in the inner part. In addition, the thin accretion disk may 
still exist in the outer parts as a so-called truncated disk emit- 
ting thermal MIR radiation. Later, observations gave indeed ev- 
idence for this accretion mode in LLAGN, and thus it is favored 
for these sources (see Ho 2008; Yuan 2007, for an overview). 



Of course a change in the accretion mode could also invoke a 
change in the UV or X-ray emission properties, which in course 
could by itself lead to an observed X-ray-deficit (or MIR-excess) 
compared to Seyferts accreting in a different mode. 

To summarize, circum-nuclear star formation, the putative 
torus, a jet or the accretion disk, can all significantly contribute 
to the observed MIR emission of LLAGN. By characterizing the 
MIR spectral energy distribution, it might be possible to distin- 
guish between these different possible components, or at least 
constrain their relative strength - and thus test the underlying 
models for LLAGN. This work represents the first step by inves- 
tigating whether the apparent deviation from the luminosity cor- 
relation of G+09 still holds for a larger sample of LLAGN, and 
finding possible explanations for any deviation. Therefore, new 
VISIR imaging data of a sample of 17 LLAGN is presented here 
and will be compared to the previous higher luminosity AGN 
sample from G+09. 

The structure of this paper is as follows: First, the sample se- 
lection and properties are described in Section 2. The following 
section (3) details the VISIR observations, MIR data reduction 
and photometric measurements. The obtained results are pre- 
sented in Sect. 4. And combined with additional objects from 
the literature the MIR-X-ray correlation for LLAGN is investi- 
gated. Comparisons to Spitzer/IRS and IRAS are presented and 
discussed in Sect. 5, as well as possible star formation contam- 
ination. In addition, a possible dependence of the MIR-X-ray 
luminosity ratio on the accretion rate is discussed in Sect. 5.4. 
The paper ends with a summary and conclusions (Sect. 6). 

2. Properties of the LLAGN sample 

The sample of LLAGN was chosen with the goal of testing the 
previously found MIR-excess/X -ray-deficit. Thus, the sample is 
not suited to detect objects significantly above the correlation 
(i.e. showing a MIR deficit). In detail, 17 LLAGN were selected 
according to the following main criteria: ;) visible from Paranal, 
ii) a nuclear X-ray detection in the 2-10 keV band with a cor- 
rected L2-iokeV ^ 10 42 erg/s, Hi) classified as LINER or Seyfert, 
and iv) an expected flux fn^m ^ lOmJy assuming that the ob- 
jects follow the MIR-X-ray correlation as found by Horst et al. 
(2008) (9 objects) or that they follow the trend of MIR-excess 
(or X-ray-deficit) seen for the 3 LINERS in the same work (9 ob- 
jects). A flux of lOmJy is roughly the lowest unresolved 12//m 
flux VLT/VISIR can detect in a narrow-band filter in a reason- 
able amount of execution time (~ 1 h) under average conditions. 
In addition, there are objects for which the selection criteria ap- 
ply but which have sufficient VISIR N-band imaging data avail- 
able in the literature (e.g. Centaurus A, M 87). These objects will 
be included in the latter part of this analysis. Thus, the total sam- 
ple comprises 26 LLAGN with 1 1 LINERs, 2 type 1 Seyferts and 
13 type 2 Seyfert. In addition, NGC 1404 which is classified as a 
narrow emission-line galaxy (NELG) was included in this sam- 
ple because it was expected to be bright enough in the MIR. In 
this work only the following main classes will be distinguished: 
type 1 (1.0-1.5), type 2 (1.8-2.0 and lh, Hi) and LINERs. The 
AGN classification of Veron-Cetty & Veron (2010) have been 
used in all cases where available (see table 2). NGC 3125, NGC 
5363, NGC 5813 and NGC 7626 are not in this catalog but have 
been classified as LINERs by Carrillo et al. (1999) and Ho et al. 
(1995). Note that some objects have multiple different classifi- 
cations (as either Seyfert or LINER) in the literature (e.g. NGC 
1097 or NGC 7213). All these classifications are defined by op- 
tical properties (emission line ratios), which can be heavily af- 
fected by extinction. Furthermore, a significant number of AGN 
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might even be completely obscured in the optical and thus re- 
main unidentified as active (e.g. Goulding & Alexander 2009). 

Additional notes on individual LLAGN included in this 
work, in particular concerning their X-ray properties, can be 
found in Appendix B. 

2.1. Distances 

As the selected LLAGN are very nearby in a cosmological sense 
(~ 20Mpc) the motion due to local gravity potentials can domi- 
nate over the Hubble flow. Thus, for most nearby galaxies it can 
be necessary to use redshift-independent distance measurement 
methods. Here, three different distance collections are used: a) 
for comparison the redshift-based luminosity distances D L from 
NED using redshifts corrected for the Earth motion relative to 
the CMB reference frame with Hq = 73, O m = 0.27, and 
f^vac = 0.73, b) the distances given by Tully (1988) employ- 
ing the Tully-Fisher (TF) relation, and c) distances of various 
methods with preference given to the most precise and recent 
measurements, with the highest priority given to direct Cepheid- 
based distance measurements and next using the surface bright- 
ness fluctuation method (see Jacoby et al. (1992) for a review and 
comparison). If not available, the values based on the TF method 
or D L were used. Note that there is no Tully-Fisher distance 
available for NGC 1667, NGC 3125 and NGC 7626. References 
for the individual sources are listed in Table. 2. However, unless 
noted otherwise distance set c) will be used. 

2.2. X-ray luminosities 

The X-ray luminosities were selected from the literature and are 
corrected for the distances adopted in each case. Preference was 
given to the most recent observations with X-ray telescopes hav- 
ing a high angular resolution, a large effective area, broad band 
energy coverage or the best combination of these. The values 
used are listed in Table 2 and represent in all cases the intrinsic 
luminosities taking into account absorption. Therefore, they can 
of course depend strongly on the spectral models used. This af- 
fects mostly Compton-thick (CT, AT H > 1.5 ■ 10 24 ctrT 2 ) objects, 
and a number of LLAGN turned out to be such CT candidates. 
Their X-ray spectra are often dominated by scattered or extended 
emission and thus appear to be unabsorbed with a much lower 
luminosity than predicted by CT obscuration. Depending on the 
model assumed this leads to intrinsic 2-10keV luminosity esti- 
mates differing by orders of magnitude (e.g. for NGC 7743). 

To identify the "true" nature of the Compton-thick candi- 
dates, it is useful to compare the estimated X-ray luminosity to 
other intrinsic indicators, such as luminosities of the forbidden 
emission lines [OIII] i5007 A and [OIV] A25.89 fim. In Panessa 
et al. (2006) a correlation between the 2-10keV and [OIII] lumi- 
nosities has been found, which can be used to predict the intrin- 
sic X-ray luminosity for the Compton-thick candidates. There 
is also evidence for significant absorption of [OIII] in the nar- 
row line region (NLR; Haas et al. 2005). Whereas, the [OIV] is 
almost unabsorbed, making it a powerful tool (Melendez et al. 
2008). Diamond-Stanic et al. (2009) present a large sample of 
local Seyferts with measured [OIV] and observed X-ray fluxes, 
whose luminosities correlate well for the 27 unabsorbed sources 
(Spearman rank p = 0.84 with log p = -7.3). A bisector fit 
(Isobe et al. 1990) to these objects yields: 

lo / Wev\ = _ Q44 ± 3 ^ + ± 00 J%l| (1) 

\ ergs 1 / \ergs l ) 



Both luminosity correlations, X-ray- [OIII] and X-ray- [OIV], are 
used to predict the intrinsic X-ray luminosities of the CT can- 
didates, and thus to distinguish between the Compton-thin or - 
thick scenario. 

In addition, it is well known that AGN can show huge X- 
ray variation, either intrinsic or due to a change of obscura- 
tion/absorption in the innermost region (e.g. Murphy et al. 2007; 
Risaliti et al. 2007). To account for this an error of factor 2 
(0.3 dex) is assumed in all cases where variation is evident or 
only one high spatial resolution X-ray observation is available, 
while 0.6 dex is applied for CT candidates. In particular, it is 
important to note that X-ray luminosities of the following ob- 
jects have been significantly revised since the sample was ini- 
tially selected: for NGC 5813 only ROSAT data were avail- 
able inferring a much higher X-ray luminosity (logL2-iokeV = 
42.1, Schwope et al. (2000)) than later found by Chandra 
(logL 2 -iokeV = 38.77), similarly for NGC 7590 ASCA data 
yielded logL 2 -iokeV = 40.78 (Bassani et al. 1999) and later 
XMM-Newton logL 2 -iokeV = 39.70 (Shu et al. 2010a). This led 
to a severe over-estimation of the MIR flux and thus an under- 
estimation of the VISIR exposure times for these objects. 

3. Observations and data reduction 

For each target a set of two filters was chosen, such that the band- 
pass of one contained only continuum emission around 12//m 
while the other contained the polycyclic aromatic hydrocarbon 
(PAH) emission feature at 1 1 .3 /im. As all objects are at very low 
redshifts «z> = 0.006), the Nellrefl (A c = 12.27 //m, half-band 
width = 0.18jum) and PAH2 (Ac = 11.25jUm, half -band width 
= 0.59 fim) filters were used in all cases. 

The VLT/VISIR observations were carried out in service 
mode between April and September 2009 (ESO period P83). 
Most observations fulfilled the required observing conditions 
(airmass < 1.3, optical seeing < 078, MIR photometric sensi- 
tivity within 20% of expected values). An exception is IC 1459, 
which was observed in conditions that resulted in a much worse 
sensitivity. Note that NGC 676 was observed again in P86 with 
longer exposure times and the more sensitive PAH2ref2 filter 
(Ac = 1 1 .88 jum, half -band width = 0.37 fim). These data will be 
used here instead of the Nellrefl data from P83. 

Standard chopping and nodding observing mode with a chop 
throw of 8" on the small field of view (19T2X 19T2; 07075/pixel) 
was performed. Therefore, all four beams fell on the detector - 
for science targets the 2 positive on top of each other (parallel 
mode), ensuring the highest possible signal-to-noise (S/N) for 
this observing mode. After consulting Spitzer/IRAC 8yum im- 
ages (or 2MASS if the former was not available) the chopping 
throw and angle for each individual object were chosen such that 
overlap of the nucleus with non-nuclear emission regions was 
avoided. The flux level of the diffuse host emission at a distance 
of 8"drops to less than 10% compared to the nucleus in 80% 
of the cases. Thus, possible contamination due to chopping in- 
duced overlaps is estimated to be < 10%. The most severe ex- 
ception is NGC 7590, which shows very complex extended MIR 
morphologies in IRAC. Thus, it was observed with the generic 
nodding template using offset of 30", but it remained undetected 
with VISIR. For each science target a photometric standard star 
was observed within 2 hours beforehand or afterwards. The ob- 
servational data are summarized in Table 1 . 

The observations and MIR properties of the archival objects 
are described in detail in the literature (Reunanen et al. (2010) 
for NGC 1386 and NGC 4486; G+09 for ESO 005-G004; Horst 
et al. (2008) for NGC 4303, NGC 4472, NGC 4698, and NGC 
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5128; and Horst et al. (2006) for NGC 4579), except NGC 1052 
(programme ID: 382.A-0604, PI: Treister). For all those objects 
the available raw data were re-reduced in exactly the same way 
as described below. 

Data reductions were performed using the ESO delivered 
pipeline. Afterwards, 2-dimensional Gaussians were fitted to all 
beams of science and calibration target under the assumption that 
all targets appear point-like. This will be justified in Sect. 4.1. 
By subtracting the constant term from the fit it is possible to ex- 
tract the unresolved flux from the local background. Thus, this 
method is not heavily affected by the background as long as it 
is locally smooth. Due to the chosen chop throw, emission struc- 
tures extended over more than 4", e.g. the diffuse emission of the 
host galaxy or weak extended star formation regions, overlap and 
are at least partially subtracted. In particular for the MIR-faint 
galaxies targeted here, detecting extended sources with VISIR 
is difficult. Additional complications arise due to fast local sky 
background variations of the atmosphere in the MIR and various 
instrument artifacts that are not understood (see VISIR manual 
Sect. 4.9). 

The absolute photometric calibration was obtained by ob- 
serving standard stars from Cohen et al. (1999) and computing 
conversion factors from the integrated intensity of the beams. 
Its precision is limited by the uncertainties in the fluxes of the 
standard stars (< 10%). In addition, a statistical error can be 
derived by comparing the Gaussian fits of the single beams for 
each target, and is given in Table 1 . Note that comparison with 
the aperture method used in Horst et al. (2008) yields consistent 
results within the uncertainties. 

For non-detected nuclei the following procedure was per- 
formed: A 2D Gaussian with the FWHM of the standard star and 
an amplitude of 2cr BG of the local background (50x50 pixel) in 
the image center was integrated resulting in the flux upper limit. 
This method was tested and verified by simulating similar beams 
on a large number of artificial and real VISIR images. In all cases 
the artificial beam was clearly visible and thus these calculated 
values represent upper limits with a confidence of > 99%. In 
addition, they are in accordance with expectations based on the 
exposure time and the sensitivity measured from the standard 
star during the observations. 

4. Results 

4. 1 . MIR morphology and photometry 

Of the 17 observed LLAGN, 7 were detected, all in both fil- 
ters. The other 10 and NGC 1404 were not detected in any filter. 
Corresponding linearly scaled images in the PAH2 and Nellref 1 
filters are presented in Fig. 1 . These images are created by co- 
adding sub-images (50 x 50 pixels) of the negative beams (with 
changed sign) to the sub-image of the central double beam. In 
general, all objects show a point-like central source with a de- 
tection significance of ~ 5cr BG in most cases. Here, <x BG signi- 
fies the standard variation of the local background (4"). The only 
exception is NGC 3125 which seems to be extended in Nellref 1 
and also in PAH2, but still the emission is acceptably fitted by 
a Gaussian. On the other hand, the detection is at a low cr-level 
(~ 3), which further complicates a reliable measurement of its 
extension. The corresponding FWHM and fluxes as well as up- 
per limits for the non-detections can be found in Table 1. For 
NGC 1667, individual beams were detected only at a 2 cr BG - 
level, but all 3 beams are located at the expected positions with 
respect to the chopping parameters. Thus, the detection signifi- 
cance rises to 3.7cr BG for the smoothed combined image (as dis- 



played in Fig. 1). However, the measured fluxes of NGC 1667 
will formally be treated as upper limits throughout this work. 
NGC 676 was not detected but a foreground star is superposed 
5" south of the galaxy center. The latter was detected in both 
filters (5cr BG , 8 mJy in PAH2ref2) constraining the flux of the 
AGN to 6 mJy in PAh2ref2 (3cr BG upper limit). For NGC 4941 
and NGC 7213 the Airy rings are evident in the PAH2 images 
indicating that the diffraction limit was reached during the expo- 
sure in those cases. 

Within the pointing accuracy of VISIR all point-source po- 
sitions are consistent with the ones derived for the galaxy nuclei 
based on other observations (e.g. 2MASS). Thus, it is inferred 
that the observed point-sources are in fact the MIR imprint of 
the LLAGN. In addition, no other emission sources of any na- 
ture are detected within 8" around the nucleus. Furthermore, the 
science targets show a ~ 20% larger FWHM compared to the 
standards. This fact is related to the much longer exposure times 
of the science observations which lead to a beam widening due to 
two different processes: MIR seeing variations on minute time- 
scales, and instabilities in the VISIR PSF due to instrumental 
effects as described in Horst et al. (2009). Both effects can lead 
to a significant variation of the FWHM during long integration 
times. In principle, this can be quantified and corrected by com- 
paring and discarding the worst chopping raw frames. However, 
in the case of the LLAGN, this is not possible, as they are too 
faint to be seen in the individual raw frames. Therefore, compar- 
isons between standard and science observations are of limited 
use even if both were performed under similar atmospheric con- 
ditions. 

In summary, there is no obvious evidence for extended emis- 
sion (except NGC 3125), non-nuclear discrete sources or elon- 
gation in any LLAGN, but all the observed MIR emission likely 
originates in a region unresolved by VLT/VISIR. 

4.2. MIR-X-ray relation 

In the following, the P83 data are combined with the archival 
data described above (Sect. 3) to include all LLAGN detected 
by VLT/VISIR so far. In order to compare the positions of the 
LLAGN in the luminosity plane with respect to the MIR-X-ray 
correlation the monochromatic 12/im continuum luminosities 
/lL^(12yum) (hereafter Lmir) are computed by using the Nellrefl 
fluxes. For this, a flat MIR SED (in AL A ) around 12/zm was as- 
sumed and no TT-correction was applied. Due to the low redshifts 
of the sample the latter effect is negligible. As already described 
in Sect. 2.2, the intrinsic hard X-ray luminosities L 2 _iok e v (here- 
after L x ) are inferred from the literature. The corresponding lu- 
minosity values can be found in Table 2 for the distance set c). 

The MIR-X-ray correlation strength is measured by using the 
Spearman Rank correlation coefficient, p, and the corresponding 
null hypothesis probability, p (r.correlate in IDL). To quan- 
tify the correlation, linear regression in logarithmic space is per- 
formed. Due to the presence of a variety of systematic uncer- 
tainties (see discussion in G+09) and a significant number of up- 
per limits it is useful to compare different fitting methods. Here, 
the Bayesian based linmix_err (Kelly 2007) and the canonical 
fitexy (Press et al. 1992), will be used. Both algorithms include 
treatment of errors on both axes. In addition, linmix_err can 
handle upper limits and intrinsic scatter as well, which makes 
it the superior algorithm for this investigation. It uses a Markov 
chain Monte Carlo method to draw random parameter sets from 
the probability distributions given by the measured data. The 
maxima of the resulting distributions of parameter draws rep- 
resent the best-fit values (details and comparison tests can be 
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Object 


Obs. date 


Filter 


1 a 
/I rest 


Exp. time* 


FWHM ["] 


Flux 




(YY-MM-DD) 




jum 


[s] 


STD C 


Obj. 


[mJy] 


IC 1459 


09-06-05 


PAH2 


11.18 ±0.59 


362 


0.29 




< 20 






Nellrefl 


12.20 ±0.18 


1443 


0.32 




< 30 


NGC 676 


10-12-22 


PAH2 


11.19 ±0.59 


1448 


0.32 




<9 






PAH2ref2 


11.82 ±0.37 


1448 


0.33 




<6 


NGC 1097 


09-07-14 


PAH2 


11.20 ±0.59 


997 


0.39 


0.49 


60 ± 10 




06-08-07 


Nellrefl 


12.22 ±0.18 


2524 


0.32 


0.40 


27 ± 4 


NGC 1404 


09-07-15 


PAH2 


11.18 ±0.59 


544 


0.57 




< 35 




09-09-06 


Nellrefl 


12.19 ±0.18 


1992 


0.38 




< 18 


NGC 1566 


09-09-06 


PAH2 


11.19 ±0.59 


362 


0.37 


0.46 


101 ± 12 






Nellrefl 


12.21 ±0.18 


1443 


0.33 


0.42 


95 ± 13 


NGC 1667 


09-09-06 


PAH2 


11.08 ±0.58 


362 


0.37 




< 16 






Nellrefl 


12.09 ±0.18 


1803 


0.38 




< 15 


NGC 3125 


09-06-17 


PAH2 


11.21 ± 0.58 


902 


0.30 


0.72 


45 ± 15 




09-06-16 


Nellrefl 


12.22 ±0.18 


1992 


0.32 


0.74 


65 ± 13 


NGC 3312 


09-07-15 


PAH2 


11.14 ±0.59 


725 


0.38 




<24 






Nellrefl 


12.15 ± 0.18 










NGC 4235 


09-06-07 


PAH2 


11.16 ±0.59 


362 


0.29 


0.38 


38 ±6 






Nellrefl 


12.17 ±0.18 


1443 


0.33 


0.37 


33 ±4 


NGC 4261 


09-06-16 


PAH2 


11.17 ± 0.59 


544 


0.30 


0.39 


12 ±2 






Nellrefl 


12.18 ±0.18 


1992 


0.36 


0.48 


15 ±2 


NGC 4594 (M 104) 


09-06-17 


PAH2 


11.21 ± 0.59 


362 


0.30 




< 11 






Nellrefl 


12.23 ±0.18 


902 


0.33 




< 17 


NGC 4941 


09-05-11 


PAH2 


11.21 ±0.59 


725 


0.36 


0.33 


68 ±5 




06-04-19 


Nellrefl 


12.22 ±0.18 


1620 


0.35 


0.35 


75 ±4 


NGC 5363 


09-05-10 


PAH2 


11.21 ± 0.59 


544 


0.33 




< 11 






Nellrefl 


12.22 ±0.18 


1443 


0.34 




< 18 


NGC 5813 


09-05-10 


PAH2 


11.18 ±0.59 


181 


0.33 




< 21 






Nellrefl 


12.19 ±0.18 


181 


0.34 




<51 


NGC 7213 


09-05-17 


PAH2 


11.18 ±0.59 


181 


0.28 


0.30 


232 ± 5 






Nellrefl 


12.20 ±0.18 


181 


0.31 


0.32 


228 ± 15 


NGC 7590 


09-06-04 


PAH2 


11.19 ±0.59 


362 


0.32 




< 16 






Nellrefl 


12.21 ±0.18 


902 


0.34 




<28 


NGC 7626 


09-06-04 


PAH2 


11.12 ±0.58 


181 


0.32 




<24 






Nellrefl 


12.13 ±0.18 


721 


0.34 




<26 


NGC 7743 


09-07-05 


PAH2 


11.19 ±0.59 


362 


0.29 




< 11 






Nellrefl 


12.20 ±0.18 


1443 


0.32 




< 13 



" : central filter wavelength in the restframe with half-band width 
* : on-source exposure time excluding any overheads 
c : FWHM of the photometric standard star 

found in Kelly 2007). Here, typically 10 5 such random draws 
will be used. Apart from an estimate of the intrinsic scatter <r m , 
linmix_err also returns a linear correlation strength, compa- 
rable to the Spearman Rank. On the other hand, the Spearman 
Rank and the fitexy will also be stated for easier comparison 
with previous work. Note that the X-ray luminosity will always 
be treated as the independent variable. This is physically moti- 
vated, as it originates in the innermost region around the black 
hole directly tracing the accretion, while the MIR emission is 
originating likely by UV heated dust further out. 



4.2.1. MIR-X-ray relation for LLAGN 

The absorption-corrected X-ray versus MIR luminosities for all 
LLAGN are displayed with filled symbols (or arrows in case 
of upper limits) in Fig. 2 together with the AGN from G+09. 



To indicate the sensitivity limits due to the observational lay- 
out (see Sect. 2), hatched areas are displayed in the plot as well. 
Objects in those areas could not be detected at a significant level 
in the P83 VISIR programme. Their borders are determined for 
each X-ray luminosity by simulating the choice of an exposure 
time for an hypothetical point-source at a representative dis- 
tance. This, in turn, results in a MIR luminosity for which the 
detection significance would drop below 3ct B g. These areas rep- 
resent only a rough guideline and are not necessarily valid for 
every individual object, especially in those cases where the ex- 
posure times were too short due to outdated X-ray luminosities 
(see Sect. 2.2). The latter explains why some upper limits (e.g. 
NGC 5813 and NGC 7626) are so far above the sensitivity limits. 

Although the LLAGN sample does not extend over a large 
luminosity range, the MIR and intrinsic X-ray luminosities are 
highly correlated. However, 2 outliers become evident by ex- 
hibiting offsets > 1 dex towards high MIR (or low X-ray) lu- 
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Table 2. Luminosities and other properties for the whole LLAGN sample. 



Object 


AGN type 


D 


Ref. 


r 


logA'H 


logL 2 _10keV 


log XL Uy m 


logM BH 


Ref. 


log /l E dd 






[Mpc] 




[pc] 


cm -2 


[erg/s] 


[erg/s] 


[M ] 






NGC 676 


Sy 2.0 


19.5 ± 8.7 


18 


33 


> 24.3 


40.79 ± 0.60 


< 40.84 


7.56 


10 


< -4.0 


NGC 1052 


L 


18.0 ± 2.4 


12 


31 


23.1 


41.18 + 0.30 


42.13 + 0.07 


8.19 


10 


-3.8 


NGC 1097 


L 


19.1 ±4.2 


19 


32 


20.4 


40.88 + 0.15 


41.23 + 0.06 


8.08 


13 


-4.2 


NGC 1386 


Sy 2.0 


18.3 ± 1.4 


14 


31 


> 24.3 


42.10 + 0.60 


42.26 ± 0.01 


7.78 


11 


-2.9 


NGC 1404 


16LG 


18.3 ± 1.4 


14 


31 


< 22.0 


39.92 ± 0.30 


< 41.25 


8.50 


11 


<-5.2 


NGC 1566 


Sy 1.5 


14.3 + 5.9 


1 


24 


<22.0 


41.62 + 0.30 


41.76 ± 0.06 


7.01 


11 


-2.6 


NGC 1667 


Sy 2.0 


62.6 ± 4.3 


16 


106 


>24.3 


42.33 ± 0.60 


< 42.25 


7.81 


11 


<-2.8 


ESO 005-G004 


Sy 2.0 


22.4 ± 10.0 


18 


38 


24.0 


41.82 + 0.30 


42.00 ± 0.08 


7.89 


20 


-3.2 


NGC 3125 


L 


19.8 ± 1.4 


16 


34 


21.7 


40.14 + 0.30 


41.87 + 0.08 


5.77 


5 


-2.0 


NGC 3312 


L 


44.7 ± 3.1 


16 


76 


< 22.0 


41.10 + 0.60 


< 42.15 


8.28 


11 


<-3.9 


NGC 4235 


Sy 1.2 


35.2 ± 15.7 


18 


60 


21.2 


41.61 ± 0.30 


42.08 ± 0.05 


7.60 


10 


-3.0 


NGC 4261 


L 


29.4 ± 2.7 


12 


50 


23.2 


41.00 + 0.30 


41.58 + 0.05 


8.72 


6 


-4.7 


NGC 4303 (M61) 


Sy2.0 


15.2 ± 3.1 


18 


26 


< 22.0 


39.04 ± 0.30 


40.72 ± 0.23 


6.52 


10 


-3.9 


NGC 4472 (M 49) 


Sy2.0 


17.1 + 0.6 


15 


29 


< 22.0 


39.00 ± 0.30 


< 41.05 


8.81 


10 


<-6.0 


NGC 4486 (M 87) 


L 


17.1 ± 0.6 


15 


29 


22.6 


40.73 ± 0.30 


41.26 + 0.05 


9.82 


9 


-6.1 


NGC 4579 (M 58) 


L 


16.8 + 3.4 


18 


29 


21.7 


41.18 + 0.15 


41.70 + 0.06 


7.77 


10 


-3.6 


NGC 4594 (M 104) 


Sy 1.9 


9.1 + 0.8 


12 


15 


21.3 


39.99 ± 0.30 


< 40.61 


8.46 


10 


<-5.4 


NGC 4698 


Sy 2.0 


23.8 ± 6.4 


17 


40 


> 24.3 


40.82 ± 0.60 


< 41.30 


7.57 


10 


<-3.8 


NGC 4941 


Sy 2.0 


19.8 ± 1.4 


16 


34 


23.7 


41.43 + 0.30 


41.94 + 0.02 


6.91 


11 


-2.5 


NGC 5 128 (Cen A) 


Sy 2.0 


3.4 ± 0.2 


7 


6 


23.0 


41.69 + 0.15 


41.71+0.01 


7.74 


4 


-3.3 


NGC 5363 


L 


19.4+ 1.4 


16 


33 


> 24.3 


41.44 + 0.60 


< 41.30 


8.36 


10 


< -4.2 


NGC 5813 


L 


32.1 ± 2.8 


2 


54 


21.1 


38.77 ± 0.30 


< 42.19 


8.53 


10 


< -5.3 


NGC 7213 


L 


21.2+ 1.5 


16 


36 


20.3 


42.09 + 0.15 


42.48 ± 0.03 


7.74 


11 


-2.7 


IC 1459 


Sy2.0 


30.3 ± 4.2 


2 


51 


22.1 


40.61 ± 0.30 


< 41.91 


9.41 


3 


<-5.4 


NGC 7590 


Sy 2.0 


27.4 ± 6.0 


19 


47 


<22.0 


39.77 ± 0.60 


< 41.79 


6.71 


8 


<-3.2 


NGC 7626 


L 


41.9 + 2.9 


16 


71 


<22.0 


39.56 ± 0.30 


< 42.13 


8.78 


10 


<-5.2 


NGC 7743 


Sy 2.0 


19.2+ 1.6 


12 


33 


> 24.3 


41.27 ± 0.60 


< 41.15 


6.64 


10 


< -2.7 



- Notes: References on distances and black hole masses. (1) average of Tully (1988), Willick et al. (1997) and D L from NED; (2) Blakeslee et al. 
(2001); (3) Cappellari et al. (2002); (4) Cappellari et al. (2009); (5) Dudik et al. (2005); (6) Ferrarese et al. (1996); (7) Ferrarese et al. (2007); (8) 
Garcia-Rissmann et al. (2005); (9) Gebhardt et al. (2011); (10) Ho et al. (2009); (11) Hyperleda database (Paturel et al. 2003); (12) Jensen et al. 
(2003); (13) Lewis & Eracleous (2006); (14) Madore et al. (1999); (15) Mei et al. (2007); (16) D L from NED; (17) Springob et al. (2007); (18) 
Tully (1988); (19) Willick et al. (1997); (29) Winter et al. (2010). r : spatial scale corresponding to the theoretical diffraction limit for a 8.2 m VLT 
telescope at 12/jm. log A = log(L Bo i/iEdd) is the Eddington ratio. 



minosities: NGC 3125 and NGC 4303 (solid grey symbols in 
Fig. 2). Both will be excluded form the correlation and re- 
gression analysis which will be justified in Sect. 5.5. For the 
luminosities of the other LLAGN, the Spearman Rank coeffi- 
cient varies between 0.78 and 0.88 (corresponding to a null- 
hypothesis probability -2.54 > log p > -3.78 ), depending on 
the distance set used. The LLAGN luminosities have been fit- 
ted with the two algorithms described above, and the fitexy 
fit to all detected LLAGN is displayed as dashed line in Fig. 2. 
The results of both fitting algorithms are in general agreement 
(compare Table 3) while linmix_err gives systematically flatter 
slopes and larger uncertainties for the parameters for these sam- 
ples. Note that the fits including upper limits are likely biased 
due to due to the outdated X-ray luminosities leading to many 
very high upper limits. Furthermore, the fitexy fit parameters 
for the different distance sets do not all agree within the calcu- 
lated uncertainties. Thus, in general, the slope of the MIR-X-ray 
correlation for LLAGN is only constrained to be between ~0.9 
and 1.1. Note that although the Compton- thickness of some ob- 
jects (e.g. NGC 7590 and NGC 7743) remains uncertain, they do 
not significantly affect the results. This is because they are not 
detected with VISIR and their upper limits do not constrain the 
results with either high (absorbed) or low (unabsorbed) intrinsic 
X-ray luminosities. 



4.2.2. Comparison to brighter AGN 

The main goal of this investigation is to determine whether 
LLAGN deviate from the MIR-X-ray correlation found for 
brighter local AGN (G+09). The latter are plotted in Fig. 2 
with empty symbols, while "well-resolved" AGN are in addition 
marked with a central black-filled circle (see G+09 for details). 
These "well-resolved" local Seyferts are probably the least con- 
taminated by non-AGN emission, and will be used for compar- 
ison with the LLAGN. Their sample properties are repeated in 
Table 3 and the corresponding power-law fit is as well displayed 
in Fig. 2 as a dotted line. In general, the LLAGN are offset by 
~ 0.3 dex from the best fit of the brighter AGN towards higher 
MIR (or lower X-ray) luminosities. In addition, the slopes of the 
LLAGN fits are all flatter (in the X-ray vs MIR plane), and in 
particular for the Di and TF distance sets. However, the uncer- 
tainties of the LLAGN fit parameters are so large in all other 
cases, that the l-cr confidence intervals overlap with the brighter 
AGN fits. Thus, there is no large systematic offset between the 
LLAGN and AGN. Furthermore, the 6 detected LINERs con- 
strain any possible difference of this class to the Seyfert class to 
be < 1 dex in the MIR-X-ray luminosity plane. 

4.2.3. MIR-X-ray correlation for all AGN 

The previous section showed that the LLAGN do not differ sig- 
nificantly from brighter AGN. Thus, it makes sense to investi- 
gate the properties of the combined AGN sample (LLAGN + all 
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Fig. 1. VISIR images of the P83 LLAGN in the filters PAH2 and 
Nellrefl. Each row shows one object sorted by right ascension. White 
black corresponds to the mean background value «BG)) and white 
color corresponds to (BG> + 10cr BG with a linear scaling. The images 
of NGC 1667 were smoothed to increase the visibility of this low-S/N 
detection. 



AGN from G+09). The corresponding correlation and regres- 
sion properties are displayed in the bottom part of Table 3. The 
correlation strength measured by the Spearman rank coefficient 
significantly increases from p = 0.88 to 0.92, and becomes in- 
sensitive to the choice of the distance set or even inclusion of the 
outliers. The same is true also for the fitexy parameters, which 
are consistent within the uncertainties for all distance sets and 
sample combinations. In particular, the fit with distance set c) is 
indistinguishable from the one for all AGN in G+09 (displayed 
as solid line in Fig. 2). Comparing the fitting algorithms fitexy 
and linmix_err, the later gives again systematically lower val- 
ues for the fitted parameters, but consistent with the former. 

The MIR-X-ray correlation is also present in the flux-flux 
space (see Fig. 3). The intrinsic 2-10 keV flux is calculated with 
the luminosities and distances given in Table 2. Here, the for- 
mal correlation strength is weaker (p = 0.72; log p = -8.0), due 
to the decreased range in orders of magnitude compared to the 
luminosity space. The slope flux correlation is very similar to 
that of the luminosity correlation for the combined sample (us- 
ing fitexy): 



log(^) = (14.35±0.82)+(1. 12+0.08) log ( ^ X(2) 

\ mJy / \erg s 'cm 1 ) 

whereas again no /f-correction is applied. The strong flux corre- 
lation and its similar slope prove that the luminosity correlation 
is not a consequence of a selection bias related to distance. 

Summarizing, there is no significant evidence for a deviation 
of the LLAGN as a group. Instead, they are consistent with the 
MIR-X-ray luminosity correlation found in previous work for 
brighter AGN (e.g. G+09). Thus, the main result of this section 
is, that the MIR-X-ray correlation extends unchanged down to 
luminosities of < 10 41 erg/s with a small scatter. 



4.2.4. Comparison to starburst galaxies 

It is interesting to compare the MIR-X-ray properties of the 
AGN to the ones of typical starburst galaxies (SB) in which the 
MIR emission is dominated by star formation. Here, the galaxy 
sample of Ranalli et al. (2003) is used, consisting of 22 typi- 
cal nearby starburst galaxies, all possessing measured 2-10keV 
X-ray luminosities presented in the same work. The hard X-ray 
emission in these objects originates mainly in X-ray binaries and 
shock heated diffuse gas. As the star formation in these galax- 
ies is occurring on global and not (only) nuclear scales, large 
aperture photometry is used for the comparison (IRAS fluxes at 
12yum). The SBs are displayed as orange stars in Fig. 2 while 
the corresponding power-law fit is denoted as the dot-dashed 
line with parameters listed in Table 3. No luminosity uncertain- 
ties are given in Ranalli et al. (2003), thus the bisector algo- 
rithm was used. The SB form a strong correlation (p = 0.92; 
log p = -9.0) with a high luminosity ratio (r = 3.1) and a small 
scatter (o> = 0.25) in agreement with Krabbe et al. (2001). They 
form a population clearly distinct from the AGN. Interestingly, 
the slope (0.97 ± 0.07) of the power-law fit is comparable to the 
one of the AGN, although the X-ray emission originates from 
completely different mechanisms. On the other hand, the nor- 
malization of the MIR-X-ray correlation for the SB galaxies dif- 
fers by ~ 2.6 orders of magnitude, towards higher MIR (or lower 
X-ray) emission. 
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Fig. 2. Absorption-corrected hard X-ray lumi- 
nosities vs. the nuclear MIR luminosities for 
all LLAGN (filled symbols) and the AGN from 
G+09 (empty symbols); blue squares: type 1 
Seyferts (type 1.5 or smaller), red diamonds: 
type 2 Seyferts; green triangles: LINERs; ma- 
genta arrow: NGC 1404 (NELG); Circinus is 
plotted for comparison only (framed in black); 
objects marked with a C: Compton thick AGNs 
(X-ray Nh > 1.5 ■ 10 24 cnT 2 ); objects high- 
lighted with central black-filled circles: "well- 
resolved" AGN from G+09; NGC 3125 and 
NGC 4303 are displayed in grey; filled or- 
ange stars: starburst galaxies (total luminosi- 
ties) from Ranalli et al. (2003); solid line: 
power-law fit to all LLAGN and AGN dis- 
played in color dashed line: power-law fit to all 
LLAGN (except NGC 3125 and NGC 4303); 
dotted line: power-law fit to "well-resolved" 
AGN from G+09; dot-dashed line: power-law 
fit to all SB. The hatched areas indicate re- 
gions below the detection limit of S/N=3 for the 
performed VISIR Nellref 1 imaging: light grey: 
objects at 20Mpc; dark grey objects at lOMpc; 
see text for details. 



5. Discussion 

5.1 . Comparison with IRAS 

The obtained high spatial resolution data give the opportunity 
to study the AGN contribution to the total 12 pm emission of 
nearby galaxies. The latter can be measured by using the large- 
aperture IRAS photometry (FWHM ~ 30"), taken either from 
Sanders et al. (2003) or the NED database. Fig. 4 shows the ratio 
of the 12 pm luminosities measured on nuclear scale with VISIR 
and global scale with IRAS over the absorption corrected X-ray 
luminosities for all LLAGN and the AGN from G+09. In gen- 
eral, the nuclear fraction of the total MIR luminosity is increas- 
ing with increasing X-ray luminosity but with large scatter of a 
factor ~ 2.5. The 2 objects showing the lowest ratio, NGC 1097 
and NGC 4303, are both dusty spiral galaxies, the former domi- 
nated by its well-known circum-nuclear starburst ring. Thus, as 
expected, the total MIR emission of the galaxies is dominated by 
non-nuclear emission and the AGN contribution dominates only 



at X-ray luminosities > 10 43 6 erg/s, while at Lx ~ 10 42 erg/s the 
host contribution to the MIR is already 3 to 10 times larger than 
the one from the nucleus. This was already found by Vasudevan 
et al. (2010) in the logarithmic space, where a strong correlation 
(p = 0.70, log p = -6.4) is present, and emphasizes once more 
the need for high spatial resolution to isolate the AGN from sur- 
rounding processes. 

5.2. Comparison to Spitzer/IRS 

Most of the LLAGN have also been observed with Spitzer/IRS, 
which enables a more detailed comparison of the innermost 4". 
In particular, low spectral resolution N-band spectra will be used 
here. Post-BCD spectra are available in the archive, created by 
a point-source profile extraction. They are well-calibrated (un- 
certainty < 10%) and are sufficient for a rough comparison with 
the nuclear properties, which are the focus of this paper. It is im- 
portant to note, however, that these IRS spectra can significantly 
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Sample 


Method 


N 




p 




logp 




0"mt 






a 




b 


r 


07 


G+09(well-resolved) 


fitexy 


22 




0.93 




-9.52 










0.19 


+ 0.05 


1.11 


±0.07 


0.15 


0.23 


G+09(well-resolved) 


linmix_err 


22 





99 ±0 


02 







12 ± 





15 


0.18 


±0.06 


1.08 


±0.09 


0.15 


0.23 


LLAGN(detected) 


fitexy 


12 




0.78 




-2.54 










0.42 


±0.27 


1.04 


±0.17 


0.40 


0.25 


LLAOJN (detected) 


linmix_err 


12 





92 ± 


16 







20 ± 





25 


0.16 


± 0.44 


0.85 


± 0.28 


0.40 


0.25 


LLAGN(all) 


linmix_err 


25 





98 ± 


05 







25 ± 





36 


0.44 


± 0.46 


1.09 


± 0.27 


0.40 


0.25 


LLAoiN (detected; D L ) 


fitexy 


12 




0.88 




-3.78 










0.06 


± 0.12 


0.74 


± 0.08 


0.42 


0.25 


LLAGNfdetected; TF) 


fitexy 


12 




0.84 




-3.14 










0.15 


± 0.17 


0.84 


± 0.10 


0.42 


0.25 


LINER(detected) 


fitexy 


6 




0.93 




-2.12 










0.49 


± 0.32 


1.01 


± 0.17 


0.55 


0.20 


LiiNriix^aeiecLeaj 


linmix_err 





u 


oy ± U 






u 


34 ± 


2 


1 




± 1 




± U.o 1 


U.J J 


u.zu 


LINER(all) 


linmix_err 


10 





94 ±0 


24 







38 ± 


5 


39 


0.39 


±2.00 


0.96 


± 1.10 


0.55 


0.20 


G+09(all) 


fitexy 


42 




0.88 




-14.01 










0.41 


±0.03 


1.12 


±0.04 


0.30 


0.36 


G+09(all) 


linmix_err 


42 





95 ±0 


03 







29 ± 





19 


0.35 


±0.06 


1.00 


±0.08 


0.30 


0.36 


LLAGN(detected) & G+09(all) 


fitexy 


48 




0.92 




-19.02 










0.41 


±0.03 


1.12 


±0.04 


0.31 


0.35 


LLAGN(detected) & G+09(all) 


linmix_err 


48 





96 ±0 


02 







26 ± 





17 


0.35 


±0.06 


1.02 


±0.06 


0.31 


0.35 


LLAGN(all) & G+09(all) 


linmix_err 


61 





98 ±0 


01 







27 ± 





17 


0.34 


±0.06 


1.06 


±0.06 


0.31 


0.35 


SB 


bisector 


22 




0.92 




-8.97 










2.96 


±0.20 


0.97 


±0.07 


3.06 


0.25 



-Notes: Methods: see text for a description; for 'linmix_err' the given values are medians of the parameters; all uncertainties are given as standard 
deviations; N: number of objects used for the analysis; p: linear correlation coefficient (for 'fitexy' : Spearman Rank); p: null-hypothesis probability; 
cr iut : intrinsic scatter (estimated by 'linmix_err'); a, b: fitting parameters of logL MIR - 43 = a + b(\ogL x - 43); F: average of the luminosity ratio 
r = log(L M i R /L x ) with o~ r its standard deviation. 
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log f,(12(im) [mjy] 



Fig. 3. Correlation of MIR and absorption-corrected hard X-ray fluxes 
for the extended LLAGN sample and AGN from G+09. Description is 
similar to Fig. 2. 



deviate from other studies, due to 2 possible reasons: in the case 
of extended objects the extraction profile under-weights the ex- 
tended emission, or the flux might be significantly lower due 
to bigger extraction areas (e.g. 20" in Gallimore et al. (2010)). 
Nevertheless, for the objects in common, NGC 1097 and NGC 
1566, the post-BCD spectra are in general agreement with those 
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Fig. 4. Fraction of total 12 pm luminosity that is emitted by the nucleus. 
It is measured by the ratio of VISIR to IRAS luminosities versus the 
absorption-corrected hard X-ray luminosity for the extended LLAGN 
sample and AGN from G+09. Description is similar to Fig. 2, except 
the dashed line, which is a power-law fit (linmix_err) to the combined 
AGN sample. 



published by Diamond-Stanic & Rieke (2010), where very small 
extraction areas have been used. 

The data are displayed in Fig. 5 while, the flux uncertainties 
of the IRS spectra (< 10%) are not displayed to maintain clarity. 
Common MIR emission lines are indicated by vertical dashed 
lines, from left to right: [Arlll] (8.99 pan), H 2 (9.67 //m), [SIV] 
(10.5 fan), PAH (1 1.3 /mi) and [Nell] (12.81 pm). For NGC 676, 
NGC 3312, and NGC 5363 no IRS spectra are available in the 
archive. Furthermore, for NGC 1667, NGC 4594, NGC 7590 and 
NGC 7626, the IRS spectra are in high resolution mode with a 
low S/N, although in some cases a strong PAH feature (1 1.3 pm) 
is still evident. In addition, archival VISIR imaging data along 
with the photometric measurements of Ramos Almeida et al. 
(2009) an (Mason et al. 2007) are shown for comparison (gray 
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data points). They show good agreement with the new measure- 
ments within the uncertainties. Only NGC 7213 shows different 
fluxes (> 10%) in both Nellrefl and PAH2 over a time of ~3 
years. The reason for this is unknown but might indicate vari- 
ability. 

For comparison with the VISIR photometry (F^ mR ), the IRS 
spectra are convolved with the normalized filter transfer func- 
tions of PAH2 and Nellrefl. The resulting fluxes, ^ RS (PAH2) 
and F^ RS (NeIIrefl), are displayed in Table 4 for all LLAGN 
detected in both filters. There are many cases, where Spitzer 
displays a much higher 12 pm continuum flux than VISIR at a 
higher spatial resolution, especially for NGC 1097, NGC 4261 
(also NGC 4486). This is also true for the upper limits derived 
for NGC 1667, NGC 4594 and NGC 7590 and implies that large 
amounts of hot dust are distributed in extended diffuse regions 
around the nucleus. Due to their faintness, the latter are proba- 
bly resolved-out in VISIR and thus not evident. Unfortunately, 
the S/N of the images is not high enough to constrain these emis- 
sion regions well, e.g., in the case of NGC 4261, the missing flux 
of ~ 35 mJy extended over a region of 0.75" diameter would 
become undetectable in VISIR. However, for the other objects 
(NGC 1566, NGC 3125, NGC 4235 and NGC 4941) the VISIR 
and IRS continuum fluxes are very similar, indicating that in 
these cases no additional significant MIR source is surrounding 
the nucleus at arcsec scale. 

The imaging in the second filter setting (PAH2) was per- 
formed to constrain the strength of the PAH emission feature 
at 1 1.3 pm, which is commonly used as a star formation indica- 
tor (e.g. Diamond-Stanic & Rieke 201 1). In particular, for NGC 
1566 the convolved IRS flux (^ RS (PAH2) = 103 + lOmJy) 
is very similar to VISIR PAH2 flux ( J F^ ISIR (PAH2) = 100.9 + 
11.5mJy). This indicates the presence of PAH emission inside 
the inner 0.4". On the other hand, the weak PAH emission fea- 
tures present in the IRS spectra of NGC 4235 and NGC 4941 
cannot be verified in the VISIR data, as the measurement un- 
certainties of are too large. For NGC 1097, the VISIR fluxes 
are both much lower than in the IRS spectrum - the latter 
showing very strong PAH emission. However, the PAH2 flux in 
VISIR is much higher than in Nellrefl implying that the nu- 
cleus of NGC 1097 also emits large amounts of PAH emission. 
Nevertheless in general, the lower continuum fluxes and absence 
of PAH emission in VISIR indicates, that this feature typical for 
star-formation usually does not originate in the very nucleus of 
LLAGN, but in the circum-nuclear surroundings (at scales larger 
than 100 pc). This was already noticed by Honig et al. (2010) for 
Seyfert galaxies. These diffuse emission regions are likely very 
extended, so that they are resolved-out on most VISIR images. 

However, to further constrain any star formation contamina- 
tion of the VISIR data, in the next section, the most extreme 
possible case - namely that all the PAH emission (seen in the 
IRS spectra) is originating in the unresolved nucleus - will be 
examined. 



5.3. Constraining the nuclear star formation contamination 



One of the main concerns of this investigation is whether the 
measured 12 pm continuum is contaminated by non-AGN emis- 
sion, e.g. circum-nuclear star formation. Such contamination 
would increasingly affect the measurements at lower luminosi- 
ties, i.e. as in the case of the LLAGN. It could in fact explain the 
MIR-excess that some LLAGN show. To better constrain the in- 
dividual star formation contamination, the LLAGN data can be 
compared to a typical star formation SED. Such template SED 
can be created by using, e.g., the IRS spectra of typical nearby 



starburst galaxies from Brandl et al. (2006). These authors col- 
lected a sample of 13 SB spectra with highest S/N. The latter 
have been made available electronically by the Spitzer/IRS Atlas 
project (Hernan-Caballero & Hatziminaoglou 2011). To create 
the SB template SED, these 13 SB are normalized by their PAH 
11.3 pm emission line flux. This approach is valid because the 
12 pm continuum scales with the PAH 1 1.3 pm flux in starbursts 
as demonstrated by a strong correlation for the whole sample of 
Brandl et al. (2006) (p = 0.98 with log p = -7.1): 



, /F v (NeIIrefl)\ 

log (— —J 



(11.82 ±0.91) 



(1.06 * 0.0 8)1 og( F(PAH :' 3 'f V (3) 
\ ergs _1 cm~ z / 



A comparable correlation was recently found by Diamond- 
Stanic & Rieke (2011). Note that this normalization results in 
a scatter of approximately 30% (0. 13dex) for the 12 pm flux 
in the individual spectra. The created SB template SED is then 
scaled for each individual LLAGN by its PAH 1 1.3 pm emission 
line flux. Next, the corresponding 12pm flux of the scaled SB 
is determined by convolution with the continuum filter function 
(Nellrefl in most cases). By this means, it is possible to estimate 
the maximum contamination by star formation in the measured 
continuum of the unresolved LLAGN nuclei in VISIR. 

The measured total PAH 11.3/im emission line fluxes 
F(PAHll3pm) are given in Table 5 for all detected LLAGN 
with evident PAH emission, while the typical measurement er- 
ror is < 10%. The scaled starburst template spectra for those 
LLAGN are plotted in orange in Fig. 5. The corresponding cal- 
culated 12 pm continuum fluxes Fy B (12pm) are indicated by or- 
ange open squares and stated in Table5 as well. Note that NGC 
1052 is not included because it was not observed in a VISIR 
filter containing only continuum, but instead in the B11.4 fil- 
ter. Although this filter includes the PAH emission feature in 
its band-pass, the latter is very weak in the IRS spectrum of 
NGC 1052. Furthermore, for the detected LINERs NGC 4261 
and NGC 4486, no significant PAH emission could be detected 
in the IRS spectra. Thus, any significant contribution of star for- 
mation to the VISIR data of NGC 1052, NGC 4261 and NGC 
4486 is very unlikely. On the other hand, the IRS spectra of NGC 
1097 and NGC 1566 show strong PAH emission and thus their 
estimated maximum star formation contribution c< B is very large 
(100% and 73% respectively). Note, however, that in the case of 
NGC 1097 and ESO 005-G004, F* B (12pm) was actually con- 
strained more by the VISIR PAH2 flux, as the latter was lower 
than the scaled flux of the SB template. Nevertheless, in the case 
of NGC 1097, F* B (12pm) still exceeds the measured 12pm con- 
tinuum flux in VISIR. However, the PAH emission feature is 
weak in all the other cases (NGC 1386, NGC 4235, NGC 4579, 
NGC 4941, NGC 5128 and NGC 7213), so that the maximal star 
formation contribution c< B can be constrained to often much less 
than 30%. In summary, a major contribution by star formation to 
the nuclear 12pm continuum (at 0.4"scale) can be excluded for 
the majority of LLAGN (9 of 12 cases). Note that exclusion of 
the other 3, possibly contaminated, LLAGN does not affect the 
MIR-X-ray correlation. 

However, it is very unlikely that all the seen PAH emission 
originates in the unresolved nucleus (as assumed in this section). 
Consequently, the real star formation contribution is probably 
much lower in all cases, as long as the PAH 11.3/Lan remains a 
valid tracer of the former. In particular, it is not verified that the 
hard AGN radiation not only destroys the PAH molecules but 
also quenches star formation. On the other hand, PAH molecules 
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Fig. 5. Comparison of VISIR photometry (red symbols: P83 data , gray symbols: re-reduced data from the archive) and Spitzer IRS spectra (black 
line). In addition, photometric measurements in the narrow-N filter for NGC 1097 (Mason et al. 2007), and the Si2 filter for NGC 1386 and NGC 
5128 and in the N filter for NGC 1386 are also plotted in gray for comparison (Ramos Almeida et al. 2009). Horizontal error bars correspond to 
the filter pass band. The IRS spectra of NGC 1 566, NGC 4594, NGC 1 667, NGC 7590, and NGC 7626 are smoothed for better visibility. Common 
emission lines are indicated by the dotted lines. For determination of the rest frame wavelengths redshifts from the NED database are used. In 
addition, scaled starburst template SED (displayed as orange line) are shown as well as the derived upper limit for the starburst contribution to the 
continuum filter measurement (orange symbols), see Sect. 5.3 for details. 
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Table 4. Comparison between VISIR and Spitzer/IRS flux densities. 



Object 


F y vlMK (PAH2) 


Ff s (PAH2) 


tf(PAH2) 


F visiK (NeIIrefl) 


F^fNellrefl) 


tf(Nellrefl) 




[mJy] 


[mJy] 




[mJy] 


[mJy] 




NGC 1097 


60.4 ± 10.2 


168.4 


2.8 


27.3 ± 4.2 


158.2 


5.8 


NGC 1566 


100.9 ±11.5 


103.2 


1.0 


95.4 ± 13 


94.0 


1.0 


NGC 3125 


45 ± 15.4 


60.2 


1.3 


64.5 ± 13.1 


59.2 


0.9 


NGC 4235 


37.7 ± 6 


41.1 


1.1 


32.8 ± 3.7 


39.6 


1.2 


NGC 4261 


11.8 ±2 


46.0 


3.9 


15.1 ±2 


51.2 


3.4 


NGC 4941 


67.6 ± 4.9 


68.0 


1.0 


75.3 ±3.9 


76.5 


1.0 


NGC 5128 


984.8 ±32.1 


1322.9 


1.3 


1488.7 ±35.1 


1767.7 


1.2 


NGC 7213 


232.2 ± 4.5 


297.5 


1.3 


227.7 ± 15 


263.2 


1.2 



- Notes: The given IRS fluxes are extracted from the BCD spectra convolved with the normalized VISIR filter transfer functions of PAH2 and 
Nellref 1 ; the PAH2 fluxes contain not only the PAH feature at 11.3 y/m but also the underlying continuum; i?(PAH2) and ^(Nellref 1 ) give the flux 
ratios (Ff s /Fj ISIR ) for both filters. 



Table 5. Estimation of maximum star formation contribution. 



Object 


logF(PAH11.3/mi) 
[erg/s/cm 2 ] 


F v SB (12/zm) 
[mJy] 


..SBfl 
L < 

% 


NGC 1097 


-11.95 


35.2" 


100 


NGC 1386 


-12.47 


32.3 


18 


NGC 1566 


-12.21 


70.4 


73 


ESO 005-G004 


-12.37 


36.7 r 


50 


NGC 3125 


-12.54 


32.9 


51 


NGC 4235 


-13.13 


8.4 


25 


NGC 4579 


-12.76 


20.2 


33 


NGC 4941 


-13.08 


9.6 


12 


NGC 5128 


-11.59 


297.4 


19 


NGC 7213 


-12.53 


33.3 


14 



" : cf signifies the maximum star formation contribution (in percent) 
to the nuclear 12/jm continuum flux (at 0.4" scale). 

* : (12jum) of NGC 1097 constrained by the VISIR PAH flux, see 
text for details. 

c : Ff (12jum) of ESO 005-G004 constrained by the VISIR PAH2 
flux, see text for details. 



might survive in the dense clumpy clouds inside the torus it- 
self for significant timescales. In addition, young massive star 
clusters emit hard radiation fields that could destroy the PAH 
molecules (see e.g. Snijders et al. 2006). At least, recent findings 
suggest that the PAH 11.3 /im remains on average a good star for- 
mation tracer also in circum-nuclear environments of AGN (e.g. 
Diamond-Stanic & Rieke 2010). 

5.4. Accretion rates 

Having shown that most LLAGN are not dominated by nuclear 
star formation in the MIR, the relation between the MIR emis- 
sion and other AGN properties can be studied, in particular the 
accretion rate. Here, the Eddington ratio will be used as a proxy 
of the accretion rate: = L Bo i/L^ with L Bo i the bolomet- 
ric luminosity and Z^dd = 1.26 • 10 38 (Mbh/A^g) erg/s being the 
Eddington luminosity. For estimating LboI the X-ray luminos- 
ity is often used (e.g., Vasudevan & Fabian 2007; Vasudevan 
et al. 2009) with L Bo i = 10L2-iokeV being a reasonable choice, 
whereas the exact value of the bolometric correction factor does 
not affect the main conclusions of this study. The black hole 
masses used in this analysis were compiled from the literature 
with preference given to the most recent determinations, and are 
listed in Table 2 and Table A. 1 for the AGN from G+09. In most 
cases the only available mass estimate comes from using the em- 
pirical M B h - cr t relation (Gebhardt et al. 2000; Ferrarese & 



Merritt 2000; Tremaine et al. 2002), and in those cases the given 
reference is for the stellar velocity ov Furthermore, the newest 
determination of the relation (Giiltekin et al. 2009) is used, for 
which the intrinsic scatter is estimated to be 0.44 dex. Here, this 
value will be adopted as general uncertainty for M SH . 

Fig. 6 shows the MIR-X-ray luminosity ratio, r = 
\og(Luml Lx), distribution versus the computed accretion rate. 
In general, the observed LLAGN have low accretion rates ( 
^Edd ^ -2.5), whereas the low-luminosity Seyferts overlap with 
the brighter AGN. On the other hand, the LINERs get clearly 
separated with lowest accretion rates. The latter is believed to be 
a critical parameter determining the AGN structure. As already 
mentioned, for low accretion rates several changes in the AGN 
structure are suggested, e.g., for the accretion disk. In addition, 
Falcke & Markoff (2000) argued that the jet might dominate the 
AGN emission at low accretion rates and presented a jet model 
that explains the observed SED of Sgr A*. A similar model was 
used by Markoff et al. (2008) to successfully fit the extensive 
multiwavelength data of M 81, a well-known LLAGN. These 
changes in the AGN structure will also be reflected in the mid- 
infrared. Although neither the accretion disk nor the jet base can 
be resolved with VLT/VISIR, if any of the two dominates the 
nuclear MIR emission, evidence should be present in the data. 
In particular, Yuan & Cui (2005) predict, that the jet (if present) 
should dominate the X-ray emission at rates below Xgt\ as 10" 5 . 
This is for example the case for NGC 4486, NGC 4594 and IC 
1459, while NGC 1052, NGC 4261 and NGC 4579 should be 
ADAF dominated. All these objects are also in the sample of 
Yuan et al. (2009), who successfully test the above prediction by 
SED fitting. In the latter work no torus model is included and in- 
stead the MIR emission is mainly produced by the ADAF. Then, 
the MIR-X-ray luminosity ratio is predicted to be roughly one 
for the ADAF dominated sources and significantly larger for the 
jet dominated sources NGC 4486 and NGC 4594. But this dif- 
ference in the ratios is not evident in Fig. 6. 

In general, the average luminosity ratio is higher for LLAGN 
(r = 0.40 + 0.25) compared to the values for the well-resolved 
(0.15 + 0.23) and for all AGN from G+09 (0.30 ± 0.36; see 
also Table 3). In particular, LINERs have a high luminosity ra- 
tio {f = 0.55 ± 0.20) compared to the type 1 (r = 0.34 ± 0.29) 
and the type 2 Seyferts (7 = 0.23 ± 0.38). However, due to the 
large scatter any differences remain statistically insignificant. As 
expected, NGC 4303 and NGC 3125 are located at peculiar po- 
sitions in this graph, especially the latter possesses an atypically 
high accretion rate for a LINER (^Edd = 10~ 2 ). In summary, 
despite the different object types the MIR-X-ray luminosity ra- 
tio can be described as constant for all AGN with a scatter of 
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Fig. 6. 12 pm to 2-10keV luminosity ratio ver- 
sus Eddington ratio logi Edd = log(L Bo i/LEdd) 
for all AGN with black hole mass estimates 
available. Symbols are similar to Fig. 2. 



0.35 dex over at least 4 (possibly 6) orders of magnitude in ac- 
cretion rate (1(T 6 < A Edd < 1). 

5.5. The outliers 

Two of the objects from the original LLAGN sample exhibit very 
peculiar positions in the different plots. They are detected only 
very weakly at ~ 3<x BG and were excluded from the correlation 
and regression analysis, because of evidence that their observed 
fluxes at multiple wavelengths are heavily dominated by non- 
AGN processes. These two objects are now discussed separately. 

5.5.1. NGC3125 

The X-ray analysis of Chandra data published by Dudik et al. 
(2005) indicated an AGN in this object, while Zhang et al. (2009) 
reports a point-like X-ray nucleus embedded in soft emission 
based on the same data. On the other hand, this galaxy is known 
to harbor young and extreme Wolf-Rayet star clusters, two of 
which have been confirmed by optical and near-infrared obser- 
vations (Hadfield & Crowther 2006), one directly its galactic nu- 
cleus. In addition, no broad optical emission lines have been de- 
tected and it is optically classified as starburst (Kewley et al. 
2001). Thus, it is unclear, if NGC 3125 really harbors an AGN, 
or if the X-ray emission originates from X-ray binaries or an 
ULX inside or near the star cluster. The MIR emission may ei- 
ther originate from dust heated by an AGN or by the central 
star cluster, on which the VISIR imaging was centered. The lat- 
ter is favored by the fact that this object is also the only spa- 
tially resolved one in the LLAGN sample. Note that the sec- 
ond and fainter star cluster at ~ 10" to the East was not de- 
tected. Comparison of the VISIR photometry and the IRS spec- 
trum shows good agreement in the continuum indicating that the 
emission seen by IRS originates from a compact nucleus. On the 
other hand, the strong PAH emission feature in the IRS spec- 
trum seems to be absent in the VISIR data. By assuming that the 
PAH 11.3yum emission is still present in the nucleus the maxi- 
mum contribution of star formation c< B is only < 51%. However, 
complex molecules are probably also destroyed in close proxim- 



ity of extremely massive young star forming regions. Evidence is 
strong that the Wolf-Rayet star cluster dominates the MIR emis- 
sion of NGC 3 125, being the reason for its peculiar position with 
respect to the luminosity correlations. Hence, the observed MIR 
flux should be regarded as an upper limit for any dusty torus 
emission of a putative AGN in NGC 3125. 



5.5.2. NGC 4303 - M 61 

This object was claimed to be undetected in Horst et al. (2008). 
A careful re-analysis of the Nellrefl image revealed a very 
weak detection at a low <x ~ 2 (compare Fig. 7). Similar to 
NGC 1667, the negative beams were found at the expected po- 
sitions. Furthermore, the detection significance increases in the 
smoothed co-added image. Still, the MIR values calculated for 
this object are unreliable due to the very low S/N. As the flux 
of this object is very low, it will unfortunately be very difficult 
to improve the S/N with any further observations. NGC 4303 
was classified as having a HII nucleus by Ho et al. (1995), al- 
though the optical line ratios are close to the AGN regime. Its X- 
ray luminosity is determined by Tzanavaris & Georgantopoulos 
(2007) based as well on Chandra data. They conclude that this 
object is a good AGN candidate even with such a low luminos- 
ity and without a detected Fe Ka emission line. Jimenez -Bailon 
et al. (2003) derive a 0.3 dex lower value from the same data. 
Both works report a hard X-ray core surrounded by a softer 
emission environment possibly associated with a circum-nuclear 
starburst ring. Colina et al. (2002) find evidence for a young star 
cluster present in the nucleus. The latter could contribute signif- 
icantly or even dominate the MIR emission at such low lumi- 
nosities. This would explain the objects extreme position in the 
luminosity plane far off the correlation. Unfortunately, there is 
no IRS spectrum covering the N-band region. As for NGC 3125 
the nuclear star cluster might easily dominate the observed MIR 
emission leading to the observed excess. Note that off-nuclear 
young clusters with comparable MIR luminosities are found for 
example in NGC 1365 (Galliano et al. 2008) and would be sim- 
ilarly unresolved at the distance of NGC 4303. 
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Fig. 7. VISIR image of NGC 4303 in Nellrefl. Left side: co-added im- 
age (50 x 50 pixel); right side: the same image but smoothed before 
co-adding; white color corresponds to (BG) + 5<r BG with a linear scal- 
ing 



6. Conclusions 

VLT/VISIR mid-infrared (N-band) imaging of 17 nearby 
LLAGN (and one NELG) is presented at a spatial resolution of 
~ 0735 in two narrow-band filters. For the 7 detections, the im- 
ages show point-like sources, and no extended emission is found 
nearby the nucleus. 

Possible sources of this unresolved MIR emission are nu- 
clear star formation, the putative dusty torus, the accretion disk, 
the jet or the narrow-line region. Any significant contribution to 
the 12 jt/m continuum by the latter is considered unlikely (Groves 
et al. 2006; Mor et al. 2009). The other possibilities are investi- 
gated by 3 different methods: 

Firstly, comparison to the MIR-X-ray luminosity correlation 
found for brighter AGN: Using X-ray luminosities taken from 
the literature, all detections plus 9 additional LLAGN from the 
archive are compared to the sample from G+09. The MIR-excess 
(or X-ray deficit) of the LLAGN is not statistically significant. 
Instead, the MIR-X-ray correlation from G+09 is also valid for 
the LLAGN, and the latter extends down to luminosities of the 
order of of < 10 41 erg/s with a small scatter. It is noteworthy that 
the local starburst galaxies with globally measured MIR and X- 
ray luminosities show a similar correlation with a comparable 
slope close to 1 but with ~ 2.6 orders of magnitude difference 
in normalization. The 2 peculiar outliers, NGC 3125 and NGC 
4303, are explained by dominant nuclear star clusters. 

Secondly, comparison to IRAS and Spitzer/IRS shows that 
the IRAS photometry is dominated by the host galaxy in most 
cases, with decreasing AGN contribution with decreasing X- 
ray luminosities. The comparison to small-aperture IRS N-band 
spectra yields much better constraints on the nuclear star for- 
mation. In particular, the PAH 11. 3 /mi emission line flux is re- 
garded as a star formation indicator. A typical starburst template 
SED scaled by the line flux observed in IRS for the individual 
LLAGN constrains the maximum 12 /mi continuum contribu- 
tion of a putative nuclear starburst. For 2 objects, NGC 1097 
and NGC 1566, strong PAH emission is indicated even at 0.4" 
scale. Star formation might dominate their nuclei in the MIR, but 
in general the weakness or absence of the PAH feature in 75% 
of the LLAGN restricts the nuclear non-AGN contamination to 
< 30%, often much less, and thus disfavors star formation as the 
dominating MIR source in LLAGN. 

Thirdly, for low accretion rates a change in the accretion 
structure of the AGN is proposed by several works, which should 
also affect the MIR properties of the LLAGN. However, the 
MIR-X-ray luminosity ratio is independent of the accretion rate, 



over at least 4 (possibly 6) orders of magnitude down to Aem ~ 
10 \ 

These results do not prove the presence of a dusty torus, but 
they are consistent with its existence in all observed LLAGN. 
Thus, the unification model might still hold for LLAGN, as al- 
ready suggested by, e.g., Panessa et al. (2006) and Maoz (2007). 
Furthermore, Elitzur & Ho (2009) recently determined a thresh- 
old for the disappearance of the torus at very low accretion rates 
depending on the bolometric luminosity. None of the detected 
LLAGN is below this threshold. On the other hand, if torus 
would be absent and the MIR emission would be dominated by 
the jet or ADAF, then there is a fortuitous agreement in both 
slope and normalization with the correlation for brighter (torus- 
dominated) Seyferts. In any case, the fact that all AGN (Seyferts 
and LINERs) can be well-described by a single correlation with- 
out any strong offsets makes the MIR-X-ray correlation a very 
useful observational tool for converting between MIR and X-ray 
powers, irrespective of the nature of the individual AGN. 

Increasing the sample size, in particular for the LINERs, 
would lead to much tighter constraints on the nature of the MIR 
emission in LLAGN. Unfortunately, current MIR facilities ei- 
ther lack the sensitivity or spatial resolution necessary to observe 
most LLAGN. On the other hand, further investigations using 
additional multiwavelength data will yield further clues on the 
nature of the MIR emission and properties of LLAGN. 
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Appendix A: Brighter AGN properties 

For completeness the properties of the AGN sample from G+09 
relevant here are stated in Table A. 1 . 
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Table A.l. Luminosities and other properties for the G+09 minus LLAGN. 
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- Notes: References for the black hole masses. (1) Bettoni et al. (2003); (2) Garcia-Rissmann et al. (2005); (3) Hyperleda database (Paturel et al. 
2003); (4) Lodato & Bertin (2003); (5) Lewis & Eracleous (2006); (6) McElroy (1995); (7) Peterson et al. (2004); (8) Wang & Zhang (2007); (9) 
Woo & Urry (2002); (10) Winter et al. (2010); (11) Woo et al. (2010). log A = log(L Bo i/L Edd ) is the Eddington ratio. 



Appendix B: Additional notes on individual objects 

8.7. IC 1459 

Detailed X-ray modeling of Chandra observations of this ob- 
ject are presented in Gonzalez-Martin et al. (2009b). Satyapal 
et al. (2004) give comparable model parameters for the same 
data set. For the X-ray luminosity the mean of both works is 
used. There is no [OIII] A5QQ1 A or [OIV] /I25.89yum luminos- 
ity measurement in the literature. The X-ray data suggest that the 
object is not Compton-thick. In addition, a compact radio core 
was detected by Slee et al. (1994), confirming the AGN nature 
of this object. Thus, the non-detection with VISIR may be ex- 
plainable by the bad weather conditions during the observations 
of IC 1459. 

B.2. NGC 676 

The only redshift-independent distance measurement for NGC 
676 available in the literature is by Tully (1988) employing the 
Tully-Fisher relation, and thus is used in this work. This ob- 
ject is included in the Palomar optical spectroscopic survey (Ho 
et al. 1995). But it was not detected in radio by Ho & Ulvestad 



(2001). The used X-ray data based on XMM-Newton is published 
in Panessa et al. (2006). They interpret this object as Compton- 
thick in accordance to Akylas & Georgantopoulos (2009) based 
on the same data. This is as well supported by the relatively 
strong [OIII] emission (Panessa et al. 2006). Unfortunately, there 
is no Spitzer data available, and thus no [OIV] measurement. 



B.3. NGC 1052 

The X-ray properties of NGC 1052 are described in Gonzalez- 
Martin et al. (2009b) based on Chandra data. The derived X- 
ray luminosity is very similar to the one given in Satyapal et al. 
(2004) for the same observation after correction for the differ- 
ent distances. The mean of both is used in this work. Counting 
as a prototypical LINER, NGC 1052 is surprisingly bright and 
obscured in X-rays. While the [OIII] luminosity is in agree- 
ment with the X-ray luminosity (Gonzalez-Martin et al. 2009a), 
there has been no [OIV] emission line detected (Satyapal et al. 
2004). Very recently Brightman & Nandra (201 1) published high 
S/N XMM-Newton data of this object stating a higher intrinsic 
logL2-iokeV = 41.5, which would bring this object closer to the 
correlation. 
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B.4. NGC 1097 

The X-ray luminosity is taken from Nemmen et al. (2006) using 
Chandra observations. It is interesting to note that the measured 
[OIV] brightness would imply a much higher X-ray luminosity 
(log Lx ~ 42.2) according to the correlation found between both 
(Diamond-Stanic et al. 2009) while the [OIII] luminosity implies 
a lower value (log L x ~ 40.0) using the [OIII]:X-ray correlation 
(Panessa et al. 2006). At larger scale (~ 10") the well-known 
starburst ring is evident in both filters but can only be partly seen 
in the small VISIR field of view. The infrared properties of the 
nucleus have already been investigated in detail by Mason et al. 
(2007). They find the absence of PAH 3.3 /im emission at subarc- 
sec scale, which seems to be typical for strong radiation fields. 
On the other hand, attempts to fit the infrared data with torus 
models failed, which led them to the conclusion that the torus 
is absent or weak. Instead, a nuclear star-forming region, co- 
existing with the LLAGN, is assumed to dominate the observed 
MIR emission. The likely PAH 11.3/im emission in the VISIR 
data indicates that the nucleus is indeed dominated by emission 
coming from such a star-forming region. This also matches with 
the finding of Storchi-Bergmann et al. (2005) of UV evidence for 
a nuclear star cluster inside the innermost 9 parsec of NGC 1097. 
But on the other hand, the MIR-excess relative to the luminos- 
ity correlation is consistent with the co-existence of star-forming 
region and a torus inside the nucleus of NGC 1097. For further 
disentangling of both components additional MIR observations 
are needed. 

B.5. NGC 1386 

NGC 1386 is one of the objects for which the different methods 
used to determine the distance show the largest discrepancies 
(~ 50 % variation). As this object is very likely Compton-thick 
it is difficult to derive a reliable X-ray luminosity. Here the ap- 
proach of Levenson et al. (2006) is followed and the luminos- 
ity is estimated by using the [OIII] line strength. The derived 
value is consistent with the Fe Ka emission line properties ob- 
served by Chandra and analyzed in the same work. It also agrees 
with the [OIV] brightness stated in Diamond-Stanic et al. (2009) 
and LaMassa et al. (2010). Thus, all available evidence favors a 
Compton-thick nature for this object. 

B.6. NGC 1404 

This elliptical galaxy is falling into the center of the Fornax clus- 
ter leading to large differences in the different distance measure- 
ment methods. It is classified as NELG in the optical. No evi- 
dence for AGN activity was found in the literature. Nevertheless, 
it is relatively bright in X-rays with a luminosity of log Lx ~ 41 
reported by Kim et al. (2006) based on Chandra observations. 
As this value is calculated for the whole galaxy including also 
soft X-ray emission it has to be regarded as an upper limit for 
any Compton-thin AGN in this object. Very recently, Grier et al. 
(201 1) re-examined the Chandra data and find a clear detection 
of a nuclear point source with log Lx ~ 40 making a LLAGN in 
this object likely. However, in this work NGC 1404 is treated as 
non-AGN galaxy and comparison object only. 

B.7. NGC 1566 

Similar to NGC 1386, NGC 1566 is showing huge discrepan- 
cies between the different distance measurement methods. The 
only published X-ray data of this objects is from ROSAT (Liu 



& Bregman 2005), where the emphasis was on off-nuclear ultra- 
luminous X-ray sources. Thus, the nuclear point-source was fit- 
ted only by a simple power-law. In the lack of any other available 
data, these values are used nevertheless. This X-ray luminosity is 
in good agreement with the predicted values from the [OIII] and 
[OIV] emission line measurements stated in Diamond-Stanic 
et al. (2009). 



B.8. NGC 1667 

Bianchi et al. (2005) finds an X-ray flux variation of a factor 
of 100 for this object and interprets the data such, that the object 
is Compton-thick. Unfortunately, they give no intrinsic X-ray lu- 
minosity or flux derived from the used XMM-Newton data. Thus, 
here the value of Panessa et al. (2006) based on ASCA observa- 
tions is used. The latter authors as well interpret the data in favor 
of Compton-thickness. This is supported by the appearance of a 
strong Fe Ka line (Bassani et al. 1999). On the other hand, the 
same ASCA data were originally used by Pappa et al. (2001) for a 
Compton-thin fit, giving a two orders of magnitude lower X-ray 
luminosity for a similar observed flux. The stated [OIII] mea- 
surements in Panessa et al. (2006) support a Compton-thick na- 
ture. The [OIV] emission (Diamond-Stanic et al. 2009; LaMassa 
et al. 2010) predicts as well an intrinsic X-ray luminosity similar 
to the one for the Compton-thick case. The VISIR observations 
support this scenario with NGC 1667 being much closer to the 
MIR-X-ray correlation than in the unabsorbed scenario. 



B.9. NGC 3312 

NGC 3312 was classified as LINER by Carrillo et al. (1999). 
The X-ray data are taken from Hudaverdi et al. (2006) which 
describes XMM-Newton observations of this object. No other X- 
ray data are available as well as no Spitzer or [OIII] observations. 



B.10. NGC 4235 

XMM-Newton data analyzed in Papadakis et al. (2008) and 
Bianchi et al. (2009) yielding similar results is used here. 
Panessa et al. (2006) derived a 0.6 dex higher luminosity but 
based on ASCA observations, which are treated with a lower pri- 
ority due to less spatial resolution. It is noteworthy that the prop- 
erties of this type 1 LLAGN are in good agreement with those of 
brighter objects, and thus it is likely that NGC 4235 accretes in 
radiatively efficient mode as well (Papadakis et al. 2008). 



B.11. NGC 4261 -3C270 

Gonzalez-Martin et al. (2009b) and Satyapal et al. (2005) arrive 
at very similar X-ray luminosities based on the same Chandra 
observation. It is interesting to note, that van der Wolk et al. 
(2010) reports a non-detection with VISIR in the SiC filter for 
this source. These observations were re-analyzed by us and a 
weak but clear detection was found. The measured flux matches 
the new observations presented here, and is 6 times higher than 
their reported upper limit (2mJy). Additional evidence for MIR 
emission of hot dust in the nucleus of this galaxy has been found 
by Leipski et al. (2009). Thus, a dusty torus is likely present in 
NGC 4261, which is in agreement with the fact that the object 
aligns well with the MIR-X-ray correlation. 
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B.12. NGC4472 

This object was included in Horst et al. (2008) erroneously 
because the nuclear X-ray luminosity was ambiguous due to 
many off-nuclear X-ray sources. Only an upper limit for the 
AGN (logLx < 39.32) derived from Chandra observations is 
given in Panessa et al. (2006). On the other hand, in a very 
detailed analysis Maccarone et al. (2003) found a nuclear X- 
ray source with logLx ~ 39 in the same data. This is in very 
good agreement with the [OIII] emission line brightness stated 
in Diamond-Stanic et al. (2009) while the [OIV] data indicate 
a much higher luminosity (logLx ~ 41.4). But the latter seems 
to be the only indication of Compton-thickness which is also not 
supported by the non-detection with VISIR. NGC 4472 might be 
a so-called "true" Seyfert 2 candidate - meaning that the broad 
emission lines are intrinsically absent rather than obscured by 
a torus (Nicastro et al. 2003). This was discussed in Akylas & 
Georgantopoulos (2009), but the XMM-Newton data presented 
there is heavily contaminated by off-nuclear emission and thus 
unreliable. Unfortunately, no conclusion about the existence of 
a torus in NGC 4472 can be drawn in this work 

B.13. NGC 4486 -M 87 

Here the mean of the X-ray luminosities of Gonzalez-Martin 
et al. (2009b) and Satyapal et al. (2004) is used. Both were us- 
ing the same Chandra data and state very comparable results. 
The flux seen in the Spitzer IRS spectrum matches well with the 
VISIR data. Interestingly, Perlman et al. (2007) explains all ob- 
served MIR emission with synchrotron emission by the jet and 
conclude that there is probably no torus present in NGC 4486. 
Very recent Herschel/PACS and SPIRE observations (Baes et al. 
2010) show the absence of significant thermal components in the 
MIR to far-infrared SED on large scale (> 10"): a single power 
law can roughly fit the overall SED from MIR all way to the 
radio with a slope of ~ -0.7. But due to the comparably low 
spatial resolution no conclusion about the presence of a torus 
can be drawn from this data. 

B.14. NGC 4594 -M 104 

X-ray properties are taken from Gonzalez-Martin et al. (2009b) 
and Ho et al. (2001) based on Chandra with marginally differ- 
ent results. This object was detected in the N-band by Grossan 
(2004) with Keck LWS. The reported nuclear flux at 10.2/im 
(11.1 mJy) is consistent with the upper limit derived here and 
would place this object close to the correlation. The Spitzer/IRS 
spectrum is very noisy and the much higher flux implies heavy 
non-nuclear contamination at this lower spatial resolution. 

Similar to NGC 4472, NGC 4594 is also is a so-called "true" 
Sy 2 candidate and was investigated in detail in Shi et al. (2010). 
They found strong silicate emission indicating that this object 
possesses a face-on dusty torus. The combined multiwavelength 
properties, including X-ray and optical/UV imply that NGC 
4594 indeed lacks broad emission lines in contradiction of the 
simple unification models. 

B.15. NGC 4698 

Cappi et al. (2006) and Gonzalez-Martin et al. (2009b) give 
X-ray spectral fits with low absorption differing by 0.4 dex de- 
rived from XMM-Newton and Chandra data (log L x = 39.39 and 
38.97), respectively. On the other hand, Gonzalez-Martin et al. 
(2009a) interpret NGC 4698 as Compton-thick using the same 



Chandra data. This conclusion is based on the [OIII] to X-ray 
flux ratio. Similarly high [OIII] fluxes are stated in Panessa et al. 
(2006) and Diamond-Stanic et al. (2009). In addition, a high X- 
ray obscuration is consistent with the [OIV] emission line bright- 
ness of the latter work. Thus, the object is very likely Compton- 
thick but with a large luminosity uncertainty. However, MIR up- 
per limit of NGC 4698 is consistent with both the Compton-thin 
or thick scenario. 

B.16. NGC 5363 

Initial X-ray studies, Gonzalez -Martin et al. (2009b) and Sansom 
et al. (2006), applied unobscured spectral fits to the available 
XMM-Newton observations and calculate a low X-ray luminos- 
ity (logLx = 39.68 and 40.30) for NGC 5363. However, in the 
most recent analysis (Gonzalez-Martin et al. 2009a) give a much 
higher luminosity from the same data but with a Compton-thick 
model. This is supported by the [OIII] brightness and a relatively 
flat spectrum above 2keV. Unfortunately, there is no Spitzer/IRS 
spectrum and thus no [OIV] measurement available. Here the 
Compton-thick fit will be adopted. 

B.17. NGC 5813 

Here the X-ray luminosity of Gonzalez-Martin et al. (2009b) 
based on Chandra is used. The X-ray image shows extremely 
diffuse soft emission. But no hard emission core is evident, 
which contradicts the existence of an AGN in this object. On 
the other hand, it could be an Compton-thick object as argued 
in Gonzalez-Martin et al. (2009a). But the [OIII] luminosity 
given in the same work predicts a only slightly higher X-ray 
luminosity than in the Compton-thin fit. The AGN nature is 
supported by the detection of a compact radio core by Nagar 
et al. (2005). But even by adopting the Compton-thick luminos- 
ity (log L x = 40.55), the non-detection with VISIR is not sur- 
prising and the derived upper limit does not constrain the corre- 
lation. 

B.18. NGC 7213 

X-ray measurements by XMM-Newton (e.g. Starling et al. 2005), 
Chandra (Shu et al. 2010b) and Suzaku (Lobban et al. 2010) 
available in the literature give very similar X-ray luminosities. 
Here, the mean value is used. As these come from different satel- 
lites and epochs, it seems unlikely that this object varies strongly 
in the X-rays. Thus, a smaller uncertainty is adopted for this ob- 
ject (0.1 dex). 

B.19. NGC 7590 

XMM-Newton data reveals that the X-ray emission of NGC 7590 
is dominated by an off-nuclear source and extended soft emis- 
sion from the host galaxy (Shu et al. 2010a). Both dominates the 
older ASCA observation by Bassani et al. (1999) with less spa- 
tial resolution. Unfortunately, for this object there is no Chandra 
observation available yet. Here the estimated nuclear hard X-ray 
luminosity of Shu et al. (2010a) is used. But note that based on 
the [OIII] to 2-10keV flux ratio the object is likely Compton- 
thick, which would imply a much higher intrinsic value. This 
is in contradiction to the hypothesis, that this object is a "true" 
Sy 2 candidate (compare also Shi et al. 2010). The upper MIR 
flux limit derived from VISIR observation does not help to con- 
strain its nature as even a 2 orders of magnitude higher intrin- 
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sic 2-10keV luminosity would be consistent with the correlation 
found. Similarly, the IRS spectrum is too noisy for further anal- 
ysis and it is unclear if the nucleus is detected here at all. 

B.20. NGC 7626 

The Chandra observation of this object has only been published 
by Humphrey & Buote (2006) but for the whole galaxy in the 
0.1-10keV band. The [OIII] flux measured by Ho et al. (1997) 
predicts a much lower value (log L x = 39.56) and will be used 
in this study. Dudik et al. (2009) analyzed the forbidden emis- 
sion lines seen by Spitzer and concluded that this object might 
not be an AGN based on the non-detection of [NeV] and broad 
Ha emission. But the upper limit for the predicted X-ray lumi- 
nosity (logLx < 41.03) from the upper-limit of the [OIV] non- 
detection is consistent with the [OIII] prediction. In addition, a 
point source was detected in radio (Nagar et al. 2005). So the nu- 
clear nature of NGC 7626 remains uncertain. The non-detection 
with VISIR favors a small X-ray luminosity but the upper-limit 
would be consistent with the MIR-X-ray correlation in any case. 

B.21. NGC 7743 

The properties of this object are very similar to NGC 5363: it 
might be either highly absorbed or not at all: based on XMM- 
Newton Gonzalez-Martin et al. (2009b) give an unobscured fit 
(log L x = 39.50), while Panessa et al. (2006), Gonzalez-Martin 
et al. (2009a), and Akylas & Georgantopoulos (2009) inter- 
pret the X-ray data as Compton-thick. The [OIII] data (e.g. 
Gonzalez-Martin et al. 2009a; Panessa et al. 2006) is favor- 
ing Compton-thickness in agreement with [OIV] measurements 
stated in Diamond-Stanic et al. (2009). Although there is no 
broad Ha line detected (Terashima et al. 2000), a compact radio 
core is present (Ho & Ulvestad 2001). In summary, the true na- 
ture of this source is unclear but likely Compton-thick. Thus, the 
X-ray luminosity given by Gonzalez-Martin et al. (2009b) will 
be used here. Note however, that a Compton-thin nature would 
be in better agreement with the MIR-X-ray correlation. 



